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ABSTRACT 


A conditional  lower  bound  on  the  minimand  of  an  integer  program  is 
a number  which  would  be  a valid  lower  bound  if  the  constraint  set  were 
amended  by  certain  inequalities,  also  called  conditional.  If  such  a 
conditional  lower  bound  exceeds  some  known  upper  bound,  then  every 
solution  better  than  the  one  corresponding  to  the  upper  bound  violates  at 
least  one  of  the  conditional  inequalities.  This  yields  a valid  disjunction, 
which  can  be  used  to  partition  the  feasible  set,  or  to  derive  a family  of 
valid  cutting  planes.  In  the  case  of  a set  covering  problem,  these  cutting 
planes  are  themselves  of  the  set  covering  type.  The  family  of  valid 
inequalities  derived  from  conditional  bounds  subsumes  as  a special  case 
the  Bellmore-Ratlif f inequalities  generated  via  involutory  bases,  but  is 
richer  than  the  latter  class  and  contains  considerably  stronger  members,  where 
strength  is  measured  by  the  number  of  positive  coefficients.  The  paper 
discusses  two  algorithms  based  on  cutting  planes  from  conditional  bounds. 

None  of  them  uses  the  simplex  method  (though  a variant  based  on  the  latter 
is  also  feasible).  Some  computational  experience  is  presented. 
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SET  COVERING  WITH  CUTTING  PLANES 
FROM  CONDITIONAL  BOUNDS 

by 

Egon  Balas 
1.  Introduction 

We  consider  the  set  covering  problem 
(SC)  min  {cx  | Ax  > e,  x^  = 0 or  1,  J e N} 

where  A = (a^)  is  mX  n»  e eRm,  e = (1,...,I),  c e Rn,  and  e {0,1}, 
ieM  = {l,...,m},  j eN  = {l,...,n}.  We  will  denote  by  a1  and  a^  the  1-th 
row  and  J-th  column  of  A,  respectively.  Without  loss  of  generality,  we 
assume  that  Cj  >0,  Vj  e N.  Using  established  terminology,  we  call  a 
vector  x satisfying  the  constraints  of  (SC)  a cover,  and  the  set  of  Indices 
j such  that  Xj  = 1,  the  support  of  the  cover.  A cover  Is  called  prime  If 
no  proper  subset  of  Its  support  defines  a cover. 

This  problem,  and  its  equality-constrained  counterpart,  the  set 
partitioning  problem,  are  useful  mathematical  models  for  a great  variety 
of  scheduling  and  other  important  real  world  problems,  like  crew  scheduling, 
truck  delivery,  tanker  routing,  information  retrieval,  fault  detection, 
stock  cutting,  offshore  drilling  platform  location,  etc.,  and  a literature 
of  considerable  size  exists  on  solution  methods  for  these  models  (see  [6] 
for  a survey  of  set  covering  and  set  partitioning;  [5]  for  a computational 
study  and  comparison  of  several  solution  techniques;  [31,  [8]  and  [9]  for 
some  of  these  methods;  and  [2]  for  a more  recent  survey  of  set  partitioning, 
which  also  contains  a bibliography  of  applications  of  both  models) . 
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In  this  paper  we  propose  a new  approach  to  set  covering,  based  on  the 
idea  of  conditional  bounds.  In  section  2 wa  introduce  this  concept  for 

arbitrary  mixed  integer  programs,  and  show  how  it  can  be  used  to  derive 
valid  disjunctions.  The  latter  in  turn  can  be  used  either  to  partition 
the  feasible  set  in  the  framework  of  a branch  and  bound  approach,  or  to 
derive  a family  of  valid  cutting  planes.  In  the  case  of  a set  covering 
problem,  the  cutting  planes  derived  from  conditional  bounds  are  themselves 
of  the  set  covering  type.  These  cuts  are  discussed  in  section  3,  where 
the  Bellmore-Ratliff  inequalities  [3]  generated  via  involutory  bases  are 
shown  to  be  a special  case  of  the  larger  family  of  cutting  planes  defined 
in  this  paper.  In  section  4 we  examine  the  conditions  under  which  a cutting 
plane  derived  from  a conditional  bound  cuts  off  a specified  prime  cover. 

Since  the  family  of  cuts  from  conditional  bounds  is  too  large  to  be  generated 
in  its  entirety,  in  section  5 we  discuss  a procedure  for  generating  "strong" 
members  of  the  family.  Sections  6 and  7 discuss  heuristics  for  generating 
"good"  prime  covers  and  feasible  solutions  to  the  dual  of  the  linear  program 
associated  with  (SC),  which  are  needed  to  generate  convenient  cuts.  Next 
we  state  two  algorithms  based  on  cutting  planes  from  conditional  bounds 
(section  8).  Section  9 contains  a numerical  example,  and  in  section  10 
we  discuss  some  early  computational  experience. 

2.  Conditional  Bounds 

The  central  idea  of  our  approach  is  to  derive  valid  inequalities  for 
the  set  covering  problem  from  conditional  bounds.  Since  this  concept  is 
meaningful  for  arbitrary  mixed  integer  programs,  we  will  introduce  it  in 
this  more  general  context. 
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A conditional  lower  bound  on  the  (objective  function)  value  of  a 
mixed  integer  program  (P)  is  a number  which  is  a valid  lower  bound  if  (P) 
is  amended  by  some  inequalities.  The  inequalities  used  to  derive  the 
conditional  bound  are  called  conditional.  The  purpose  of  these  inequalities 
is  to  produce  a conditional  lower  bound  at  least  equal  to  a known  upper 
bound.  If  this  is  achieved,  then  at  least  one  of  the  conditional  inequalities 
is  violated  by  any  solution  better  than  the  one  associated  with  the  upper 
bound,  and  this  yields  a valid  constraint. 

To  be  more  specific,  consider  the  integer  program 

(P)  min  {cx  | Ax  > b,  x > 0,  x^  integer,  J e N} , 

where  A is  an  arbitrary  mXn  matrix  (n  * |n|)  and  b is  an  arbitrary 
m-vector.  The  pair  of  dual  linear  programs  associated  with  (P)  is 


(I*) 

min  {cx  | Ax  > b. 

x > 0} 

and 

(D) 

max  {ub  | uA  < c. 

u > 0} 

For  any  feasible  optimization  problem  (S) , let  z(S)  be  the  value  of 
(an  optimal  solution  to)  (S) . 

Any  feasible  solution  x to  (P)  provides  an  upper  bound  cx  on  z(P), 
and  any  feasible  solution  u to  (D)  provides  a lower  bound  ub  on  z(D)  = z(L), 
hence  also  orf'z(P). 


Let  u be  a feasible  solution  to  (D) , i.e.,  such  that 
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(1)  uA  < c,  u > 0, 

and  suppose  the  constraints  of  (L)  and  (P)  are  amended  with  the  set  of 
(conditional)  inequalities 

(2)  Cx  > e , 

- P 

where  C * (c.  .)  is  a pXn  matrix  of  0's  and  l's,  and  e is  the  p-vector 

P 

of  l's,  with  1 < p < n.  Suppose  also  that  C has  no  zero  i^ows,  i,e., 
n 

(3)  Ec  >1  , 1»  1.....P. 

j=l 

Further,  denote  by  (P  ) and  (L  ) the  problems  obtained  by  amending 

L*  C 

the  constraint  sets  of  (P)  and  (L) , respectively,  with  the  set  (2)  of 

conditional  inequalities,  and  let  (D  ) be  the  dual  of  the  linear  program  (Lr) . 

L*  L 

Let  u be  a m-vector  satisfying  (1).  If  there  exists  a p-vector  v,  v > 0, 
v ^ 0,  such  that 

(4)  vC  < c - uA, 

then  (u,v)  is  a feasible  solution  to  (D„)  and  therefore  ub  + ve  is  a lower 

C p 

bound  on  2(0^)  = z(Lc),  hence  also  on  z(Pc>.  We  will  say  in  this  case  that 

ub  + ve  is  a conditional  lower  bound  on  z(P). 

P 

Now  let  Zy  be  a known  upper  bound  on  z(P) . If 

(5)  vep  > Zy  - ub, 

j 


/x  ^'*4* 
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i.e.,  If  the  conditional  lower  hound  on  z(P)  exceeds  or  equals  the  upper 
bound  then 

*<PC>  > ^(Lc) 

> ub  + ve  > z 

P ~ U’ 

and  hence  every  feasible  solution  to  (P)  better  than  the  one  associated 
with  the  bound  z^  violates  at  least  one  of  the  inequalities  (2),  i.e., 
satisfies  the  disjunction 

P n 

V ( Z ex  t I)- 
i=l  j=l  J J 

Since  C is  a 0-1  matrix  and  x is  a 0-1  vector,  this  disjunction  is  the 
s ame  a s 

P " 

v < E c x = 0). 

i«l  J=1  J J 
If  we  denote 

Qt  = {.  j £N  I ci  j = 1}  , i =l,...,p, 

the  result  that  we  have  just  proved  can  be  stated  as  follows. 

Theorem  2.1.  Let  z^  be  a known  upper  bound  on  z(P),  and  let  u satisfy 
(1)  . If  there  exists  a pXn  matrix  C = (c^),  c_e[0,l},  V-  i , j (1  < p < n) , 
satisfying  (3),  and  a p-vector  v > 0,  v ^ 0,  satisfying  (4)  and  (5), 
then  every  feasible  solution  x to  (P)  such  that  cx  * z^T  satisfies  the 
disjunction 
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P 

(6)  V (x  = 0,  jeQ  ). 
i*=l  J 

We  have  stated  Theorem  2.1  for  the  pure  integer  case  in  order  to 
simplify  the  exposition,  but  it  (and  the  rest  of  this  section)  is  easily 
seen  to  carry  over  to  the  case  of  a mixed  integer  program.  Indeed,  if 
N^CN  is  the  set  of  integer-constrained  variables,  4 N,  ~et  c ^ = 0, 

¥ i = l,...,p,  V-  jeN\Nj,  anc  the  results  are  valid  with  some  changes  in 
the  notation. 

Note  that  the  disjunction  (6),  though  derived  from  a conditional  bound, 
is  an  "unconcitionally"  valid  constraint. 

The  first  question  that  arises  in  connection  with  ifieoretn  2.1,  is  that  of  rh 
existence  of  a matrix  C and  a vector  v satisfying  the  above  requirements. 

Theorem  2.2.  Let  be  a known  upper  bound  on  z(P) , let  u be  a 
feasible  solution  to  (D),  and  let 

(7)  s = c - uA. 

Then  there  exists  a oair  C,  v satisfying  the  requirements  of  Theorem  2.1 

if  and  only  if 

(8)  2 s > z - ub. 

Proof.  Let  u and  s satisfy  (1)  and  (7).  If  (8)  holc^,  then  the  pair 
C,  v,  defined  by  C = I (the  identity  matrix  of  order  N)  and  v.  = s ^ , jeN, 
satisfies  (3),  (4)  and  (5). 

Conversely,  if  C,  and  v satisfy  (3),  (4)  and  (b),  then  adding  the 

inequalities  (4)  and  substituting  s for  c - uA  yields 


r » > £ v ( z c ) 

JtN  J l-l  jeN  J 

> ve  [ ircnc  (3)  ] 

“ P 

> - ub  [from  (5)]  Q.E.D. 

Note  that,  if  p * 1,  i.e.,  if  C has  a singie  row,  (hen  the  disjunction 

(6)  becomes  x.  = 0,  ieQ.  . Somewhat  more  generally,  we  iv  cl.  iollowir.;.  . 

J * 

Remark  2,1.  Let  z , u and  s be  as  in  Theorem  2.2,  an;  ce  ine 
Q0  = (jeNjs.  > zu  - ut>3- 

Then  every  feasible  solution  x to  (P)  such  that  c>:  ; tisfies 

J 

x.  = 0 , jeQQ. 

Thus,  whenever  Q t 0,  che  variables  indexed  bv  Qn  ca . oe  set  to 
0 o 

zero  permanently. 

Example  1.  Let  (?)  be  an  integer  program  (minimi  tion  f o 
10  variables,  and  let  = 35  be  the  value  oi  some  k.r.  i it  -,  r solution 
Further,  let  j be  a feasible  solution  to  the  di  ..  of  tl 
program,  with  ub  = 27.9,  an;  let  the  reduced  coses  £ . ^ sc 


wi  t'n  ui  bo 
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We  have  - ub  = 35  - 27 . 9 = 7 . 1 . To  construct  a pair  C,  v satisfying 

the  conditions  of  Theorem  2,1,  we  first  choose  v.  - s v_  - s,  and  v„  * sQ; 

1 1 L 4 jo 

then  v^  + v,  t » 7,2  ^ J.l.  Next  we  choose  the  l's  among  the  elements 

c . , of  C in  such  a wav  that  v,c, . + v c„ . + v.c_,  < s,  for  * l,. ...10. 
ij  ' 1 lj  2 2j  3 3j  — j 

(Neither  v,  nor  C is  of  course  unique.)  C and  v are  shown  oelow,  where 
the  blanks  are  0's. 


From  Theorem  2.1,  every  x satisfying  the  constraints  of  (?)  and  such 
that  cx  < 35  satisfies  the  disjunction 

X1  = X2  = x10  = 0 V x4  = x6  = x7  = x9  = 0 V X6  = xd  ’ XI0  = 0 • 

Given  a feasible  solution  u to  i'D1)  whose  associated  reduced  costs  s 
satisfy  (8),  it  is  usually  not  difficult  to  find  a pair  C,v  satisfying 
the  conditions  of  Theorem  2.1;  but  the  problem  is  to  find  one  which  defines 
a conveniently  strong  disjunction  (6).  The  following  procedure  provides  a 
general  framework  for  doing  this,  which  allows  for  manv  variants. 


1. 


Choose  a minimum-cardinality  subset  S of  N,  such  that 


*■- ■+?  'ZsszMbanK:'*-:- 
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and  order  S = {j  (1) , . . . , j (p ) } according  to  decreasing  values  of  s 
Then  set  v±  = i = 1,...  p,  and  c^.,  = 1 for  h = l , chJ(i)  = 0 

for  h t i,  i = 1 , . . . , p . 

2.  For  j e N\S,  define  recursively  for  i = l,...,p, 

f i-1 

; 0 or  1 if  E v c + v.  < s . 

| h=l  h hJ  1 ~ J 


otherwise. 


The  option  of  setting  = 1 or  = 0 in  step  2 represents  a 
choice  between  including  i :.nto  Q. , or  leaving  it  available  for  inclusion 
into  one  or  more  sets  Q This  can  be  decided  by  efficiency  criteria, 

as  will  be  seen  later.  It  easy  to  check  that  any  pair  C,  v constructed 
by  this  procedure  satisfies  the  requirements  of  Theorem  2.1. 

A disjunction  of  the  form  (6)  obtained  from  a conditional  bound  can 
be  used  to  partition  the  feasible  set  in  the  framework  of  a branch  and 
bound  procedure,  by  creating  p subprobleias  defined  by  the  constraints 


Xj  = 0 , jtQj. 


2 x > 1 ; x = 0 , j eQ, 

jcQi  J 1 


E x > 1 , 

jeQi  J 


J<Vi 


h 2 1 


For  certain  classes  of  integer  programs,  this  way  of  branching  seems 
highly  efficient.  It  is  currently  being  tested,  for  instance,  in  a new 


penalty  method  for  solving  traveling  salesman  problems  [ij,  with  excellent 


computational  results. 
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Another  way  of  using  the  disjunction  (6)  is  to  generate  from  it  a 
family  of  valid  inequalities.  In  the  case  of  a set  covering  problem, 
these  inequalities  turn  out  to  be  of  the  set  covering  type,  as  shown  in 
the  next  section. 

In  a broader  context,  the  idea  of  deriving  a valid  ("unconditional") 
constraint  from  one  or  several  conditional  constraints  may  have  many  other 
applications.  One  of  them  appears  in  [7],  where  a properly  chosen  inequality 
is  used  to  derive  a bound  from  the  fact  that  either  the  inequality  or  its 
complement  must  be  satisfiec  by  any  feasible  solution. 


3.  Cutting  Planes  From  Conditional  Bounds 


From  now  on,  we  address  ourselves  to  the  set  covering  problem  (SC) 
introduced  in  section  1;  i.t.,  A = (aij)  with  a_e{o,l},  v i,j,  and  b = e, 
where  e is  the  m-vector  of  I's.  We  will  denote 


= {j  cN  | a = 1],  j eM. 


Theorem  3.1.  Suppose  the  conditions  of  Theorem  2.1  are  satisfied, 
i.e.,  the  disjunction  (b)  is  a valid  constraint  for  (SC).  With  each 
ie[l,...,p),  associate  an  index  h(i)eM,  such  that  ^ t.  Then 

denoting 


every  cover  x such  that  cx  < satisfies 
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v • 


K 
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Proof.  For  i = 1,...,?,  the  i-th  term  of  the  disjunction  (6)  implies 


E x > 1 , Vh  eM  . 

J eN^\Q i 

Hence,  for  any  choice  of  indices  h(i)  eM,  i = 1,...,P,  (6)  implies  the 
disjunction 


V ( £ x > 1)  , 

‘-1  J‘Vi>VQi 

which  in  turn  implies  (10). 


Q.E.D. 


The  cutting  planes  of  Theorem  3.1  are  set  covering  inequalities, 
valid  in  the  sense  of  being  satisfied  by  every  cover  better  than  a given 
one.  Since  these  properties  are  the  same  as  those  of  the  Eellmore-Ratlif f 
cuts  [3]  obtained  by  the  u9e  of  involutory  bases,  we  next  examine  the 
relationship  between  the  latter  and  our  inequalities  from  conditional  bounds. 
First,  we  show  in  the  next  theorem  that  the  Bel lmore-Ratlif f inequalities 
are  a subclass  of  the  class  of  inequalities  defined  by  Theorem  3.1.  Then 
we  show  by  way  of  example  that  the  subclass  in  question  is  a proper  one. 

Theorem  3.2.  The  Bellraore-Ratliff  inequalities  [3]  are  a subclass 


of  the  class  (10) . 

Proof . Let  x be  a prime  cover,  B an  involutory  basis  associated 

with  x,  and  c - c a.  the  j-th  reduced  cost,  where  c is  the  m-vector 

j B J B 

whose  i-th  component  is  c^  if  the  basic  variable  associated  with 

row  i is  (the  structural  variable)  x^.^,  and  0 if  the  basic  variable 
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associated  with  row  i is  a slack.  (When  B is  an  involutory  basis,  the 
reduced  costs  are  known  to  be  of  the  above  form.)  The  Bellmore-Rat li f f 
cut  associated  with  x and  B is  then 


E x > 1 
JeF  J 


where 


F = [J  eN  | c - y < 0} 


To  obtain  this  cut  from  a conditional  bound,  let  1^  = {j  (1) , . . . , j (p)} 
be  the  index  set  of  the  structural  basic  variables,  and  set  u * 0,  s = c, 
and  v^  =*  j(i)  el^.  Then  u satisfies  (1)  and  v satisfies  (5)  (with 

equality)  for  z^  = cx. 

We  now  construct  the  matrix  C = (c^j)  Theorem  2.1  as  follows. 

Let  h (i ) be  the  row  index  associated  with  basic  variable  x Define 


ah(i),j  ’ jeN'J 


, j eF 


i = 1, • • • »P  > 


Qi  = W 


\F  , i = 1 , . . . ,p  • 


C trivially  satisfies  (3).  To  see  that  together  with  v it  alsc 
satisfies  (4),  note  that  for  jeN\F, 

ill  ViCij  ” ih  cjdA(«,j 


< ci  - ; 
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Thus  the  Bellmpre-Ratliff  cuts  are  a special  case  of  the  cuts  (10). 
Furthermore,  they  are  a proper  subclass  of  the  class  (10),  i.e.,  the 
conditional  bound  approach  yields  other  inequalities  besides  those  ob- 
tainable via  involutory  bases.  Some  of  those  other  inequalities  are 
considerably  stronger,  in  the  sense  of  having  fewer  positive  coefficients. 
This  is  illustrated  by  the  next  example. 

Example  2.  Consider  the  set  covering  problem  whose  costs  c^  and 
coefficient  matrix  A are  shown  in  Tableau  1. 
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1 

1 

Tableau  1. 

The  0-1  vector  x whose  support  is  {2,3,5,12,13,17}  is  a cover, 
satisfying  with  equality  all  the  inequalities  except  for  1 and  8,  which 
are  oversatisfied.  To  apply  the  Bellmore-Ratlif f procedure,  one  associates 
with  x an  involutory  basis.  The  variables  X£,  x^,  x^  can  be  basic  only  in 
rows  3,  4 and  6 respectively.  Since  rows  1 and  8 are  slack,  the  variables 
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and  x^  can  be  basic  only  in  rows  11  and  10  respectively.  Finally,  x^ 
can  be  basic  in  any  of  the  4 rows  2,  5,  7,  9;  and  accordingly  there  are  4 
involutory  bases  that  can  be  associated  with  x.  We  will  denote  them  by 
1*2,  and  B^,  according  as  x^  is  basic  in  row  2,  5,  7 or  9 

respectively.  The  basis  B2  (after  row  permutations)  is  shown  in  Tableau  2. 


3 

4 
6 

11 

10 

2 

5 

7 
9 
1 

8 


2 3 5 12  13  17  25  27  29  21  28 


Tableau  2. 


The  variables  X2^,...,x^  are  slacks.  The  basis  B,.  can  be  obtained 
from  B2  by  interchanging  rows  5 and  2,  and  replacing  the  slack  variable  X2^ 
in  row  5,  by  X22  In  row  2.  The  other  two  bases  can  be  obtained  in  an 
analogous  way.  The  4 cutting  planes  that  can  be  obtained  by  the  Bellmore- 
Ratliff  procedure,  depending  on  which  basis  is  used,  are 


X1  + X6  + x9  + x10  + x15  + X16  + X13  + x20  2 1 


x4  + x6  + xg  + x1Q  + xn  + x19 

x6  + x7  + x1Q  + x15  + x19  + x2Q 

X6  + X8  + X10  + X14  + X18  + X20 


> 1 
> 1 
> 1 


from  B„ 


from  B? 


from  B. 


from  Bg. 
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On  the  other  hand,  using  the  conditional  bound  approach,  we  construct 
(by  Inspection  or  a heuristic)  the  dual  vector 

u = (0,1, 1,1, 1,1, 2, 0,1, 2, 2) 

which,  together  with  the  associated  reduced  cost  vect''- 


s - (2, 0,0, 2, 0,0, 0,1, 1,0, 1,1, 1,1, 1,0, 0,2, 0,1) 

satisfies  the  condition  (1). 

The  cover  x whose  support  is  {2,3,5,12,13,17}  yields  z^  = cx  * 14; 

and  the  dual  vector  u yields  the  lower  bound  ue  « 12. 

Since  zy  - ue  =2,  QQ  = {j  e N | Sj  > 2}  = {1,4,18},  and  thus  (Remark  2.1) 

every  cover  better  than  x satisfies  x = x.  = x,  = 0.  Hence  we  replace  N 

I io 

by  N\{l,4,18j.  Further,  to  apply  Theorems  2.1  and  3.1,  we  set  p = 2,  with 

V1  * S12  8nd  v2  c s13’ 

812  + 913  “ 2 

> Zy  - ue  . 

Then,  using  (for  instance)  the  conditional  inequalities  (and  the  matrix  C) 


defined  by 


Ql  = {8,12,14,20}  , Q2  = {9,11,13,15} 

we  obtain  the  disjunction 

X8  = X12  = x14  = X20  = ° V x9  = X11  = X13  = X1 5 = 0 ' 
Applying  Theorem  3.1  with  h(l)  = 11,  h(2)  = 10,  we  have  (using  the 
newly  defined  set  N"> 


W«i  * (6>19) 


\(2)\02  ’ l10>16’19! 
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and  thus 

W - {6,10,16,19}  . 

We  have  obtained  the  cut 


x6  + x10  + x16  + x19  ^ 1 


which  has  only  4 positive  coefficients,  whereas  each  of  the  involutory 
basis  cuts  has  at  least  6. 

The  above  inequality  cuts  off  x.  This  is  due  to  the  way  we  chose 
the  components  of  v and  the  row  indices  h(i)  eM,  as  will  be  shown  in  the 
next  section.  If  we  drop  the  requirement  for  a specified  cover  to  be  cut 
off,  we  can  obtain  inequalities  which  are  "stronger"  in  the  sense  of 
having  fewer  positive  coefficients.  By  a Judicious  choice  of  the  column 
indices  j(i)  for  which  we  set  = s^^,  ant*  row  indices  h(i),  one 
can  generate  the  cuts  with  the  fewest  possible  positive  coefficients.  Thus, 
for  instance,  the  choice  v^  = s^»  ^ = s^,  and  h ( 1 ) = 8,  h(2)  = 5,  yields 
the  valid  inequality 

X17  + X19  ^ 1 5 

whereas  v^  = s^,  v2  •-=  s^  and  h(l)  = 8,  h(2)  = 4 yields 

x3  ^ 2?  i • 


4.  Some  Properties  of  Cuts  from  Conditional  Bounds 


In  order  to  obtain  a conditional  bound,  and  thus  to  be  able  to 
generate  the  type  of  constraints  discussed  in  the  previous  section,  one 
needs  a feasible  solution  u to  (D) , whose  associated  reduced  costs  satisfy 
(8).  This  requirement  is  easy  to  meet.  The  next  theorem  and  its  Corollary 
give  a broad  sufficient  condition  for  u to  satisfy  them. 
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Let  S(x)  denote  the  support  of  x. 

Theorem  4.1.  Let  x be  a cover  for  (SC)  and  let  u and  a aatisfy 
(1)  and  (7) . If  u also  satisfies 

u(Ax  - 1)  = 0, 

then 

(12)  2 s *=  cx  - ue  . 

JeS (x)  J 

Proof.  Let 

S+  = {j  eS(x)  | Sj  > 0} 
and  consider  the  pair  of  dual  linear  programs 


(Lx) 

min  {cx  | Ax  > e 

* XJ 

> 1,  j eS+,  Xji  o,  j e N\S+} 

and 

CDX) 

max  {ue  + 2 

JeS 

8a 

J 

! U8j  + Sj  “ cjt  j e N;  u > 0, 

s > 0}  . 

Clearly,  x is  a feasible 

solution  to  (L^),  and  (u,s)  i 

s a feasible 

solution 

to  (D^).  Further,  x 

and  (u,s)  satisfy  the  complementary 

slackness 

conditions 

(a1*  - l)ut  = 

0 

, i eM 

(Xj  - l)8j  = 

0 

, j eS+ 

xjSj  = 

0 

, j e N\S+ 

Hence  x and  (u,s)  are  optimal  solutions  to  (L^ ) and  (D.),  respectively, 


which  implies  the  equality  of  the  two  objective  function  values. 


Q.E.D, 
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1 


y 

N 

y 


•i 
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An  Immediate  consequence  of  Theorem  4.1  is  the  following. 

Let 

T(x)  = [icM|aix  = l]  . 

Corollary  4.1.1.  Let  x be  a cover,  and  let  u and  s satisfy  (1)  and 
(7).  If  u also  satisfies 

(13)  u - 0 , V-ieM\T(x), 

then  there  exists  [j  (1) , . . . , j (p) ) c S+  such  that  v defined  by  = s 
i “ l,...,p,  together  with  u,  satisfies  (5). 

Proof.  If  u satisfies  (13)  , then  together  with  x it  satisfies  the 
complementarity  condition  of  Theorem  4.1,  and  therefore  (12)  holds. 

From  (12),  v^  as  defined  in  the  Corollary,  together  with  u,  satisfies  (5). 

Q.E.D. 

Thus,  any  feasible  solution  to  (D)  that  has  positive  components  only 
for  ieT(x)  produces  a vector  v which  satisfies  (5) . A procedure  of  the 
type  outlined  in  section  2 can  be  used  to  find  a matrix  C which  satisfies 
the  other  requirements  of  Theorem  2.1.  Choosing  appropriate  rows  h(i)eM 
will  then  produce  a cut  of  the  family  (10)  defined  in  Theorem  3.1. 

From  the  point  of  view  of  having  a finitely  convergent  procedure, 
it  is  essential  that  the  inequalities  generated  be  not  only  valid  but  also 
new,  i.e.,  no  cut  should  be  generated  twice. 

Given  two  inequalities  of  the  form  (10)  with  associated  index  sets 
and  W^,  the  first  inequality  implies  the  second  one  if  a".d  only  if 

We  saY  that  an  inequality  (10)  is  new,  if  it  is  not  implied 

by  any  inequality  of  the  current  problem  (SC). 


»*<r 
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Theorem  4.2.  The  Inequality 

(10)  E x > 1 

jeW  J 


* 


of  Theorem  3.1  is  new  if  and  only  if  N\W  is  the  support  of  a cover. 

Proof . If  N\W  is  the  support  of  a cover,  say  x,  then  x violates 
(10),  whereas  it  satisfies  all  the  inequalities  of  the  current  problem  (SC). 
Hence  (10)  is  new.  Conversely,  if  N\W  does  not  define  a cover,  then  there 
exists  a row  ieM  such  that  N1  0 (N\W)  = 0,  i.e.,  N^W.  But  then  the  cut 
(10)  defined  by  W is  either  identical  to,  or  implied  by,  the  i-th  in- 
equality of  (SC).  Q-E«D. 

Theorem  4.2  gives  a necessary  and  sufficient  condition  for  an  in- 
equality (10)  to  cut  off  at  least  one  cover.  Next  we  give  a sufficient 
condition  for  an  inequality  (10)  to  cut  off  a specified  prime  cover  x. 

As  before  we  will  denote 

Mj  - {i  e M | a^  -1}  , J cN  . 

Remark  4.1.  If  x is  e prime  cover  for  (SC),  then 
1^  n TOO  t 0 , Vj  eSOO  . 

Proof.  Follows  from  the  definition  of  T(x)  and  that  of  a cover  prime. 

Theorem  4.3.  Let  x be  a prime  cover  for  (SC),  let  z^  be  a known  upper 
bound  on  the  value  of  (SC) , and  let  u and  s satisfy  (1)  and  (7) . Further, 
let 

(9)  E s a z{J  - u e 

j e S J 


1 
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(H) 


vi  = 


• 1 " 1 » • • • i P * 


and  let  C be  any  pXn  0-1  matrix  satisfying  conditions  (3),  (4)  and 

1 h = 1 

1 ~ !»•••»?  • 

0 h 4 1 


(15) 


Finally,  let  = {jeNlc^  = l}  and 


where 


<16>  h(i)eI(x)nMj(1))l=l,..„P< 

Then  the  Inequality 


(10) 


E x.  i 1 


JeW 


J 


is  satisfied  by  every  cover  x such  that  cx  < z^,  and  cuts  off  x. 

Proof.  By  assumption,  the  pair  u,  s satisfies  (1)  and  (7),  while 
the  matrix  C satisfies  (3)  and,  together  with  the  vector  v defined  by 
(14),  satisfies  (4).  From  (9),  u and  v satisfy  (5).  The  existence  of 


such  a pair  C,  v follows  from  (9)  and  Theorem  2.2.  The  additional  require- 
ment  (15)  on  C can  always  be  met.  Hence  the  pair  C,  v defined  in  the 
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r 


On  the  other  hand,  from  (15),  j(l)eQ^,  i =*  Hence 

soo  n [Nh  = 0,  i = 1,...,P> 

i.e.  S(x)f)W  = 0 and  the  inequality  (10)  cuts  off  x. 

Q.E.D. 

The  results  of  the  last  two  sections  can  be  used  to  generate  a 

family  of  inequalities  whose  members  are  different  from  each  otfier  and 

from  the  initial  constraints  of  (SC)  , and  are  satisfied  by  every  cover  whose 

value  is  smaller  than  z^.  Since  the  family  of  such  cuts  is  of  considerable 

size,  one  would  like  to  be  able  to  generate  some  of  the  "stronger"  members, 

according  to  some  reasonable  measure  of  strength.  This  will  be  discussed 

in  the  next  section.  First,  however,  we  examine  another  question  related 

to  the  properties  discussed  above. 

If  the  set  Qq  defined  in  Remark  2.1  is  nonempty,  then  the  variables 

indexed  by  Q can  be  set  to  0.  Since  the  vector  u (and  the  associated  s) 

0 

used  in  the  definition  of  Q^,  is  not  required  to  have  any  other  property 

but  to  satisfy  (1),  the  question  arises  as  to  which  vector  u is  more 

likely  to  produce  a nonempty  set  Q^,  or  more  generally,  a larger  set 

While  this  question  is  hard  to  answer  in  general,  there  is  a clear  answer 

12  12 

for  the  case  when  the  choice  is  between  two  vectors  u , u such  that  u < u . 

Theorem  h.4.  For  k = 1,2,  let  u satisfy  (1),  and  define 

V k 

s / = c - Z u a , jeN, 

J J i€M  lJ 

and 

Qg  * [ J sn| s^  > Zy  - uke}. 

12  12 

Then  u < u implies  £ Q . 


•'X  *-■ 


Proof.  If  J€Qq>  then 


1 1 
s . = c - E u a 

J J 1 lJ 


> ZU  - u e. 


c + E u > z ; 
J ieM\Mj  1 U 


and  if  u u , the  last  inequality  implies 


Ci  + Ui  - 2U  ’ 


i.e. , 


2 _ 2 
s = c - E u a.  . 

J J ieM  1 1J 


> ZU  ‘ u e, 


hence  jeQg. 


Q.E.D. 


Thus,  the  set  defined  with  respect  to  a vector  u satisfying  (1) 
is  always  contained  in  the  set  Oq  defined  with  respect  to  any  vector  u 
obtained  from  u by  increasing  some  of  its  components  (while  maintaining 


uA  < c) . Therefore,  in  generating  the  set  Q^,  one  should  always  use  a 
"maximal"  dual  vector  u,  i.e.,  one  whose  components  cannot  be  increased 
without  decreasing  some  component  or  violating  uA  < c. 


5.  Generating  Cuts 


In  this  section  we  discuss  two  procedures  for  generating  conveniently 
strong  members  of  the  family  of  cuts  introduced  above.  Since  all 


inequal ities  have  a right  hand  side  of  1 and  left  hand  side  coefficients 
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of  0 or  1,  we  will  use  as  a measure  of  strength  the  number  of  positive 
cocf f icient:.  (the  smaller  the  number,  the  stronger  the  Inequality).  'Ihe 
fact  that  ar.  inequality  A is  stronger  than  B in  this  sense,  i.e.,  has 
fewer  positive  coefficients,  does  not  imply  of  course  that  A dominates  B 
or  implies  B. 

Procedure  CUT  generates  a member  of  the  subfamily  characterized  in 

Theorem  4.3,  i.e.,  an  inequality  which  is  satisfied  by  every  solution  x 

such  that  cx  < z , and  which  also  cuts  off  a specified  prime  cover  x. 

The  strengt  . of  an  inequality  (10),  i.e.,  the  size  of  the  set  W,  depends 

on  the  into  er  p and  the  size  of  the  sets  ti . \Q. , 1 = l,...,p,  of 

nfi)  i 

Theorem  3.1.  To  have  p con/eniently  small,  the  procedure  set6  v^  = 

for  a aeque.ee  of  indices  j(i),  t = l,...,p,  corresponding  to  the  largest 

reduced  costs  s , jeS(x).  each  row  index  h(i),  i = l,...,p,  is  of  course 

chosen  from  T(x)  TiM  ^ , as  prescribed  by  Theorem  4.3.  Further,  in  order 

to  have  W as  small  as  possible,  tne  sequence  of  row  indices  a(i)  is  chosen 

so  as  to  come  as  close  as  possible  at  each  step  ke[l,...,p,  to  minimizing 

the  set  W.  \W.  , , where  WA  = 0 and 

k k*i  o 


• k 


S^nce  j S . — 1 implies  (Remurk  2.1)  that  some  variable  car.  be  permanently 
set  to  zero,  we  assume  that  this  has  already  been  done  for  all  variables 

indexed  by  Q„  , .nd  thus  there  are  no  singleton  sets  S satisfying  (9). 

J 
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Let  M and  N be  the  row  and  column  index  sets  of  the  current  problem  (SC), 
let  x be  a prime  cover  for  (SC),  and  let  S(x)  = [jeNjx^  = lj , 

T(x)=  {ieM|aix  = l}.  Further,  let  u and  s satisfy  (1),  (7)  and  (9)  for 
S = S+,  where 

S+  = {j£S(x)|Sj  > 0}. 


CUT 


Step  0.  Initialize  WQ  = 0,  Kx  = S , yx  = ue,  s = s.  Let  t - 1 

(iteration  counter)  and  go  to  1; 

Step  1.  Define 


vt  - min{zb,  - yt,  max  s*}  , K*  = {jeK^s*;  = vt} 


j £iC 


p.  = [jeNjs^  2 » M(t)  = U M 


jeK* 


Choose  i(t)  such  that 


iNi.'t>  X(Pt  U Wt-1)  • = -Win  lNiX(PtUWt-l)' 

icT(x)  flM(t) 


(breaking  ties  arbitrarily),  and  let 


- 1 


[j(t)}  = Ni(t)  OK-  . 


Then  define 


= Wt.1UlNi(t)\Pt],  yt+1  = vt  + vt. 

If  y . > z„,  go  to  2.  Otherwise,  let 
J t+ 1 — U 


S!  ’Tt 


otherwise , 


K = Kt\(jvt)},  set  t - t + 1,  and  go  to  1. 


-•* - 
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Step  2.  Define  W = and  add  to  the  constraints  of  (SC)  the  inequality 


Z x > 1 . 
jeW  ^ 

Theorer  5.1.  In  js+|  iterations,  procedure  CUT  generates  an  inequality 
satisfied  by  every  cover  x such  that  cx  < z^,  and  violated  by  x. 

Proof.  Since  K.  *=  S+  decreases  by  one  element  at  each  iteration,  the 
procedure  ends  in  js  j iterations.  Next  we  show  that  the  inequality  generated 
by  the  procedure  satisfies  the  conditions  of  Theorem  4.3. 

The  vectors  u,  s satisfy  (1)  and  (7).  The  vector  v useu  in  the  procedure 
satisfies  condition  (14)  of  Theorem  4.3,  while  the  matrix  C used  in  the 
procedure  (not  stated  explicitly)  is  defined  by 

1 ' ptnNi(t> 

0 otherwise 


This  satisfies  requirement  i31  and  also  (15),  since  {j(t)]  = N^^OK* 
is  unique  [irom  i(t)cT(x)]  end  j(t)  t j(t-l),  t = 2,...,p.  From  the  definition 
of  the  sets  P_,  it  follows  that  C and  v satisfy  (4).  The  stopping  rule 

- z(j  guarantees  that  the  set  S = (j  (1)  , . . . , j (p) } , where  p is  the  numb., 
of  iterations,  satisfies  (9;.  Finally,  the  choice  of  i(t)  in  step  1 insures 
that  condition  (16)  is  met.  Q.E.D. 

Note  that  the  above  inequality  was  obtained  for  a given  value  of  z^. 

If  later  the  upper  bound  z^.  is  improved,  i.e.,  replaced  by  z^  < z^,  the  cut 
can  be  strengthened.  Namely,  we  have 


Corol 1 ry  5 .1.1.  Suppose 


(17) 


> 1 


is  a cut  obtained  in  p iterations  of  CUT  for  an  upper  bound  Zy,  and 
subsequently  zy  is  replaced  by  z^  < z Then  every  cover  x for  which 
cx  < z^  satisfies 


(18) 


1, 


where  k is  the  greatest  integer,  if  it  exists, 
(19) 


tc 

2 v 


t:xv i-i  - zu ' v 


such  that 


or  else  k = 0. 

Proof.  At  the  end  of  p iterations  of  CUT,  we  have 


(20) 


P+1 


= ue  4-  £ v 

i=l 


^ V 


which  implies  that  the  inequality  (17)  is  satisfied  by  every  cover  better 
tnan  the  one  associated  with  z^.  Subtracting  (191  from  (20)  yields 


p-k 

Vk+1  = ue  + > z’ , 

J.—  i. 


which  implies  that  the  inequality  (18)  is  satisfied  by  every  cover  better 
than  the  or.e  associated  with  z^.  Q.E.D. 

The  information  required  in  order  to  be  able  to  strengthen  a cut 
generated  by  CUT  1 whenever  z T?  is  improved,  is  the  sequence  of  pairs 
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6.  Generating  Prime  Covers 


Any  algorithm  based  on  the  cutting  planes  introduced  in  this  paper 
needs  to  generate  repeatedly  prime  covers  for  the  current  problem  (SC), 
and  feasible  solutions  to  the  dual  (D)  of  the  linear  program  associated 
with  (SC).  In  this  section  we  discuss  some  heuristics  for  generating 
reasonably  good  prime  covers  for  (SC)  and  feasible  solutions  to  (D) , "good" 
in  the  sense  of  giving  tight  bounds. 

The  procedure  PRIMAL  1,  to  be  described  below,  is  a kind  of  "greedy" 
algorithm,  which  selects  at  each  step  the  locally  most  promising  column, 
and  in  one  pass  produces  a ^usually  remarkably  good)  prime  cover  x. 

As  before,  let  M and  N be  the  row  and  column  index  set  of  the  current 
problem  (SC),  with  = fieM|a  = 1}  and  = {_  j eN  I j = lj  • Further,  let 

z„  be  ar  upper  bound  and  z a lower  bound  on  the  value  of  (SC). 

U L 

PRIMAL  1 

Step  0.  Initialize  R^  = m,  S^  = 0,  t •-  1,  and  go  to  L. 

Step  1.  If  Rt  = 0,  set  S = St  and  go  to  2.  Otherwise,  define 


I(t)  = {i  eR  | IN  | = min  J N | } N 

hcRt 


«*>  ' uno"1 


If  N = 0.  stop;  the  cover  corresponding  to  z^  is  optimal, 
Otherwise,  choose  J(t)  such  that 

Cj(t)/|Rt°Mj(t)!  = TOln  lcj/|Rt:n*jU 

J€NI(t) 

if  this  minimum  is  uniip  . If  it  is  not,  and  it  is  attained  for 
jsj(t),  |J(t)j  > 1,  choose  j(t)  such  that 
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lRt  'Miml  = max  K nMi 

JeJ(t)  C J 


Arid  if  there  are  further  ties,  break  them  by  maximizing  j i . 

Then  set  S[t1  = St  U fj(t)},  Rt+1  = Rt\Mj(t)»  t «-  t + 1,  anc  go  to  1. 

Step  2.  Consider  the  elements  i e S in  order  of  decreasing  costs  c^, 
and  if 

■J  M = M, 
jeS\fi}  J 

remove  i from  S.  If  all  ieS  have  been  considered,  x defined  by  x^  = 1, 

j eS,  x =0  otherwise,  is  a prime  cover.  If  cx  < z , set  z = cx  and  stop. 

U u 

If  z(.  < z^,  x is  optimal. 

The  logic  behind  the  choice  of  j(t)  in  Step  1 is  that,  the  smaller 
the  number  of  columns  that  can  be  used  to  cover  a particular  row,  say  1, 
the  higher  Che  price  to  be  paid  for  not  covering  optimally  row  i is  likely 
to  be.  Therefore  the  procedure  first  chooses  a row  i(t)  with  a minimum 
number  of  l's  in  it,  then  among  the  columns  which  can  be  used  to  cover  row 
i(t),  it  selects  one  with  minimum  cost  per  row  covered.  Ties  are  broken 
by  choosing  the  column  which  covers  the  largest  number  of  new  rows  (at 
the  same  unit  cost),  or,  if  there  still  is  a tie,  the  column  which  covers 
the  largest  total  number  of  rows.  When  all  the  rows  have  been  covered 
(X  = 0)  , step  2 checks  whether  the  cover  is  prime,  and  makes  it  prime  if 
necessary  by  dropping  redundant  columns  in  order  of  decreasing  cost. 


The  above  procedure  generates  a prime  cover  for  the  current  problem 
(SC)  by  starting  from  scratch.  Though  this  procedure  is  relatively 
cheap,  it  is  nevertheless  often  desirable  to  have  a procedure  which  starts 
with  a partial  cover  rather  than  from  scratch.  For  instance,  if 
S(x)  H Qq  + 0,  where  S(x)  is  the  support  of  the  current  prime  cover  x and 

is  the  set  of  Remark  2.1,  i.e.,  if  some  variables  which  are  at  1 in  the 

current  cover  can  be  fixed  at  0,  it  is  desirable  to  find  a new  cover 

starting  with  S = S(x)\Qq.  Also,  when  x is  a prime  cover  for  the  current 

problem  (SC)  but  ceases  to  te  a cover  because  of  the  addition  to  (SC)  of 

new  cuts,  x can  be  used  as  a partial  cover  to  start  the  procedure  for  finding  a 

new  cover.  PRIMAL  2 is  a modified  version  of  PRIMAL  1,  which  differs  from 

the  latter  only  in  the  initialization  step,  where  S^  = 0 is  replaced  by 

S,  = S(x),  che  support  of  the  partial  cover  that  we  wish  to  start  with. 

Finally,  it  often  happens  that  a good  dual  solution  u is  found,  i.e., 
one  which  yields  a high  lower  bound  ue,  but  which,  together  with  the  associated 
reduced  cost  vector  s,  does  not  satisfy  (9)  for  S = S(x),  where  x is  the 
current  cover.  If  (9)  cannot  be  satisfied  for  any  set  S £ IV,  then  of  course 
we  have  no  choice  other  than  to  modify  u and  s,  i.e.,  decrease  some  coraponen. 
of  u,  which  will  then  increase  (usually  by  more)  the  sum  of  the  reduced  cost3. 
If,  however,  u and  s satisfy  (9)  for  S = N,  then  it  is  worth  trying  to  change 
the  cover  rather  than  the  vector  u;  for  if  a cover  (any  cover)  x can  be 

found  such  chat  (9)  holds  for  S = S(x),  then  a cut  can  be  generated  whose 

strength  depends  only  on  u,  s and  z^,  but  not  on  cx.  In  other  words,  no 


-31- 


matter  how  bad  the  new  cover  x is  In  term6  of  the  associated  objective 
function,  if  It  makes  (9)  hold  for  u,  s then  It  provides  for  a better 
starting  point  for  a cut  than  could  be  obtained  by  keeping  the  old  cover 
x and  changing  u and  s. 

The  procedure  PRIMAL  3 is  meant  for  situations  like  this.  It  starts 


with  the  old  cover  x,  and  successively  introduces  variables  x.  such  that 

J 

8 >0.  Every  time  such  a variable  is  introduced,  another  variable  x^ , 

such  that  s,  = 0,  is  dropped.  Furthermore,  x.  is  chosen  so  as  to  cover 
r‘  J 

at  least  one  row  not  covered  by  other  variables  with  positive  reduced 
cost,  and  x^  is  chosen  so  as  to  leave  at  least  one  row  covered  only  by  Xj . 
These  choice  rules  are  intended  to  keep  the  cover  prime  if  possible.  When 
this  procedure  cannot  be  continued  further,  then  the  partial  cover  at  hand 
is  completed  by  using  as  few  columns  as  possible,  again  in  an  attempt  to 
keep  the  cover  prime.  The  cover  obtained  this  way  is  then  made  prime, 
if  necessary,  by  removing  tome  variable(s)  x^ , preferably  such  that  s.  = 0. 
As  berore,  let  T(x)  = [ieM|a*x  = l},  S+  *=  [jeS(x)|6j  >0},  and  denote 
= {jcN|s. ..  > 0}.  We  assume  that 


L ,s , > z„  - ue. 


,r+  J - U 


PRIMAL  3 


Step  0.  Initialize  S = S(x),  K..  = N \S  , R * M\  b M , and  go  to  1. 

1 1 JeS+  J 

Step  I.  If  Rt  = 0,  set  S = St  and  go  to  3.  If  f 3 but 

M HR  = 0,  ¥ jeK  , or  s > 0,  ¥ jeS  , set  R = M\  E M.  t..d  go  to  2. 

J J t t jeSt  J 

Otherwise  choose  j(t,l)eKt  and  j(t,2)eSt\N+  such  that 
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lMj(t,l)nRtl  " “*  lMjnRtl 

J 6K-t 


M.  , (lM  . ..OR  I = max  |M.  ClM..  ..  fl  R I. 

' j(t,2)  j(t,l)  t1  J£S  XN+  ' j J(t,l)  t1 

The"  St+1  = StU{j(t,l)}\{j(t,2)},  Kt+1  = Kt\{j(t,l)},  Rt+,  = Rt\MJ(t>1), 
t *-  t + 1,  and  go  to  1. 

Step  2.  If  = 0,  set  S = St  and  go  to  3.  If  R£  / 0 but  M H = 0, 

V-  jeN\St>  stop;  the  cover  corresponding  to  is  optimal..  Otherwise,  choose 
j(t)  such  that 

iMi(t'i  nRtl  = ciax  |m  n R | , 

JV  ; j€N\Sc  J 


set  St+1  » StUCj(t)},  Rt+1  = Rt\M  t - t + 1,  and  go  to  2. 

Step  j.  Consider  first  the  elements  ieS  such  that  s^  = 0,  then  the 

remaining  ieS,  and  if 

U M = M, 

j eS\[i}  J 

remove  i from  S.  If  all  ieS  have  been  considered,  x definec  by  x . = 1, 

jeS,  Xj  = 0 otherwise,  is  a prime  cover.  If 

Z s > z - ue, 

*eS  J u 

then  s can  be  used  to  generate  a cut.  Otherwise  u and  s have  to  be  modified, 
Next  we  turn  to  the  problem  of  generating  "good"  feasible  solutions  u 


to  (D). 
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7 . Generating  Feasible  Solutions  to  (D) 

As  already  mentioned,  the  quality  of  a solution  u to  (D)  is  given  by 

the  quality  of  the  lower  bound  ue  on  the  value  of  (SC),  provided  that  u and 

the  associated  reduced  costs  s.  satisfy  (9)  for  S = S(x).  Among  two 

feasible  solutions  to  (D)  with  the  same  value  ue,  the  one  with  a higher 

value  of  o 3 I s.  is  likely  to  produce  a disjunction  (6)  with  fewer  terms 
jeN  J 

and  more  elements  in  each  term,  i.e.,  a stronger  cut.  Hence  the  purpose 
of  our  heuristics  will  be  to  approximate  as  much  as  possible  the  lexicographic 
maximum  of  the  two-component  vector  (ue,  a). 

This  is  what  the  procedure  DUAL  1,  to  be  described  below,  tries  to 
accomplish.  It  successively  assigns  values  to  components  of  u,  which  are 
maximal  subject  to  the  dual  constraints  and  the  earlier  value-assignments. 

The  order  in  which  the  values  are  assigned  is  crucial  to  the  quality  of  the 
resulting  solution,  and  DUAL  1 considers  the  rows  of  A in  order  of  increasing 
N^.  Once  a value  has  been  assigned  to  each  component  of  u,  the  set  of 
Remark  2.1  is  identified  and  the  corresponding  variables  are  set  to  zero. 

If  u and  the  associated  3 satisfy  (9)  for  S = S(x),  then  s can  be  used 
to  generate  a cut.  Otherwise  either  the  cover  must  be  changed  as  discussed 
in  the  previous  section  (see  PRIMAL  3),  or  the  vector  u muse  be  adjusted 
(see  DUAL  3 below). 

Let  S (x)  be  the  support  of  a prime  cover  for  the  current  problem  (SC), 

and  let  z and  z be  an  upper  and  a lower  bound,  respectively,  on  the  value 
U L 


of  (SC). 
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Step  0. 
Step  1. 


DUAL  1 

Initialize  R^  = M,  = 0,  s''  •=  c,  t — 1,  and  go  to  1. 

If  Rt  = 0,  go  to  3.  Otherwise  choose  i(t)  such  that 


i(t) 


and  set 


/ , t 

min  Sj  i = i(t) 


t+1  1 

Ui  = < 


jeN 


i(t) 


t+1 


otherwise , 


t , t t+I  j «N  , x 
sj  + Ul  - u.  i(t) 


otherwise 


Define 

F 


(t) |Si+1  = » M(t)  = U M 

jeFt  J 


set  R * R \M(t),  t - t + 1,  and  go  to  1. 

Step  2.  If  uCe  '•  z , set  z = ute,  store  s*"  in  place  of  the  vector  s 
L L 

associated  with  the  previous  bound  z , and  stop.  If  z £ z , the  cover 

L L o 

associated  with  z.  is  optimal. 

As  already  mentioned,  the  choice  of  i(t)  is  crucial  for  the  efficiency 
of  the  procedure.  Since  the  rule  of  Step  1 usually  leads  to  ties,  one  may 
want  to  think  of  secondary  criteria.  Several  such  criteria  are  currently 
being  tested.  One  which  has  been  found  to  frequently  improve  the  quality 
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of  the  solutions  ts  to  give  preference  to  ieT(x)  over  ieM\T(x)  in  choosing 
i(t)  in  Step  1. 

While  DUAL  1 is  a computational ly  cheap  one  pass  procedure,  it  starts 

from  scratch  every  time  a new  pair  of  vectors  u,  s is  to  be  found.  Next 

we  state  a procedure  (DUAL  2)  which  starts  from  a dual  vector  u and 

modifies  it  after  the  addition  of  a new  cut  and  the  finding  of  a new  cover, 

by  setting  to  1 the  variable  associated  with  the  new  inequality,  and 

adjusting  the  values  of  some  other  variables. 

Let  m be  the  index  of  the  last  inequality  added  to  '(SC)  after 

generating  the  current  dual  vector  u = (u, , . . . ,u  ,)  and  tne  associated 

1 m- 1 

reduced  cost  vector  s.  Further,  let  x be  a cover  for  the  current  (SC), 

S(x)  its  support,  and  T(x)  = {ieMja^x  = l}. 


DUAL  2 


Step  0.  Initialize 

1 i = m 

ui  = 0 i e[M\T (x) ] 

u otherwise, 

i 1 

1 1 - 

s = c - u A,  and  = {jeN|s^  < O}.  Set  t *-  1 and  go  to  la. 

g-te£.  la-  If  Kt  = 0>  sat  K°  = {jeNjsJ  = 0),  Rt  = £i«M|x±  OK®  = 0}, 
and  go  to  1.  Otherwise,  choose  i(t)  such  that 


K(t)nK;i 


max 


i«M\ (mj  jut>0 


n.  n k" 
i t 


set 


otherwise 


otherwise . 


I 


l 


8 . 
J 


Then  set  Kt+..  = [jeN|Sj+^  < 0},  t •-  t + 1,  and  go  to  la. 

Steps  1-2  are  the  same  as  in  DUAL  1. 

Finally,  we  have  to  deal  with  the  situation  when  a pair  u,  s has  to 
be  adjusted  so  as  to  satisfy  (9)  for  some  S = S(x) , where  x is  a given 
prime  cover.  The  following  procedure  is  meant  to  achieve  this. 


Step  0. 
and  go  to  1. 
Step  1. 


DUAL  3 

Initialize  u^  = u,  s^  = s,  U 
Choose  i(t)  such  that 


tleM\T(x) (u*  >0],  Set  t - 1 


1 Ni (t)  0 S (x)  I 

Then  set 


max 

ieU 


|n  ns(x) 


i = l(t) 
otherwise  , 


t+1 

s. 


jeN 


ift) 


otherwise 


If  ut+‘  and  stM  satisfy  (9)  for  S = S (x)  , stop:  s = 
used  in  CUT  1.  Otherwise,  set  U . = U fi(t)},  t - t + 1, 


c-f-1  , 

s can  be 

and  go  to  1. 
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The  heuristics  discussed  here  find  reasonably  good  dual  vectors  u at 
a low  computational  cost.  Another  possibility  is,  of  course,  to  solve  (D) 
by  the  simplex  method.  While  this  is  more  costly,  only  computational 
testing  can  show  whether  the  improvement  in  the  bound  ue  is  worth  the 
extra  computational  effort. 


8 . Algorithms 

In  this  section  we  present  two  algorithms  for  solving  set  covering 

problems  by  cutting  planes  from  conditional  bounds. 

They  use  as  ingredients  the  procedures  discussed  in  the  last  three 

sections  (for  generating  prime  covers,  feasible  dual  vectors  and  valid 

cutting  planes),  and  the  following  three  subroutines  for  updating  cuts  and 

fixing  variables  whenever  possible. 

, old  . new. 

STRENGTHEN  (to  be  used  whenever  PRIMAL  replaces  zy  by  ). 

In  every  cut  E x.  > 1,  replace  W by  W , , where  k is  the  greatest 

jeW  J P P_k 

, , P old  new 

integer  such  that  v +. . .+  v , ^ z - z 

e p p-k+1  — U U 

TEST  1 (to  be  used  after  STRENGTHEN). 

Let  s be  the  reduced  cost  vector  associated  with  an<^  let  S(x)  be 

the  support  of  the  last  cover  x.  Define 

% * > zy6W  * V’ 

set  x^  = 0,  JeQq.  and  N — N\Qq.  If  NtS(x),  the  cover  associated 

, , new  . ^ . . 

with  Zy  is  optimal. 


i 
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TEST  2 (to  be  used  after  DUAL). 

new 

Same  as  TEST  1,  except  for  the  fact  that  s,  z and  z are  replaced 

U Li 

by  s11,  z^T  and  ute  respectively,  where  u*"  is  the  last  dual  vector  obtained, 
s*1  is  the  reduced  cost  vector  associated  with  u , and  z^  is  the  current 
upper  bound. 

Algorithm  1 starts  by  initializing  the  upper  bound  z^  at  £c^  and 
the  lower  bound  z^  at  0.  A typical  iteration  consists  of  the  following 
sequence  of  steps. 

A.  Use  PRIMAL  1 or  2 to  generate  a prime  cover  x.  If  z^  is  improved, 
go  to  Cl;  otherwise  go  to  Bl. 

Bl.  Use  DUAL  1 to  generate  a dual  vector  u and  go  to  C2 . 

B2.  U e DUAL  3 to  adjust  u and  s and  go  to  D. 

Cl.  Apply  TEST  1 to  fix  variables.  If  NtS(x),  stop.  Otherwise, 

if  x is  still  a cover,  go  to  Bl;  else  go  to  A. 

C2 . Apply  TEST  2 to  fix  variables.  If  NGS(x),  stop.  Otherwise, 

if  x is  still  a cover  and  (9)  holds  for  S = S(x),  go  to  D;  if  (9)  does 

not  hold  for  S = S(x),  go  to  B2 ; and  if  x is  not  any  more  a cover,  go  to  A. 

D.  Use  CUT  to  generate  a valid  inequality  which  cuts  off  the  last 
cover.  Add  the  new  inequality  to  (SC)  and  go  to  A. 

Algorithm  1 seems  the  simplest  possible  way  to  use  the  cuts  from 
conditional  bounds.  Algorithm  2,  stated  below,  is  a more  sophisticated 
procedure  baaed  on  the  same  approach,  which  (a)  strengthens  the  cuts 
whenever  z is  improved;  (fc)  avoids  starting  from  sera  .ch  c\ery  time  a 

v 

new  cover  or  a new  dual  solution  is  generated  by  using  the  "updating" 
procedures  PRIMAL  2 and  DUAL  2,  with  a periodic  return  to  PRIMAL  1 


£K£Z*UKX  ’?-■ 
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r 


and  DUAL  1;  and  (c)  uses  PRIMAL  3 whenever  z Is  improved  and  (9)  does 

L 

not  hold  for  the  current  cover,  to  find  a prime  cover  x such  that  u and 
s satisfy  (9)  for  S = S(x).  Since  the  chances  for  finding  a better  cover 
than  the  current  best  one  are  steadily  diminishing  during  the  procedure, 
the  parameter  a defining  the  number  of  iterations  after  which  we  return 
to  PRIMAL  1 is  doubled  every  time  the  latter  procedure  is  applied;  while 
the  parameter  9 defining  the  frequency  of  return  to  DUAL  1 is  kept  constant. 

After  initializing  z^,  z^,  a and  3,  a typical  iteration  of  Algorithm  2 
consists  of  some  of  the  following  steps. 

Al.  Use  PRIMAL  1 to  generate  a prime  cover  x,  and  double  the 
value  of  a.  If  is  improved,  go  to  Cl;  otherwise  to  Bl. 

A2 . Use  PRIMAL  2 to  generate  a prime  cover  x.  If  z^  is  improved,  go 
to  Cl;  otherwise  go  to  B2,  or  (if  no  cut  was  added  after  obtaining  last 
vector  u,  or  B2  was  used  3 times  since  last  use  of  Bl)  to  Bi. 

A3.  Use  PRIMAL  3 to  generate  a prime  cover  x.  If  the  last  u and  s 
satisfying  (9)  for  S = S(x),  go  to  D;  otherwise  go  to  B3. 

Bl.  Use  DUAL  1 to  generate  a dual  vector  u and  go  tc  C2. 

B2.  Use  DUAL  2 to  generate  a dual  vector  u and  go  to  C2. 

B3.  Use  DUAL  3 to  adjust  u and  s and  go  to  D. 

Cl.  Use  STRENGTHEN  to  strengthen  the  cuts,  and  TEST  1 to  fix  variables. 

If  NCS(x),  stop.  Otherwise,  if  x is  still  a cover,  go  to  El;  else  go  to  A2 . 
C2.  Use  TEST  2 to  fix  variables.  If  N£S(x),  stop.  Otherwise,  if 

x is  still  a cover  and  (9)  holds  for  S = S(x),  go  to  D;  if  (9)  does  not 

hold  for  S = S(x)  but  holds  for  S = N,  go  to  A3;  and  if  (9)  does  not  hold 

for  S = N,  go  to  Bl;  finally,  if  x is  not  any  more  a cover,  go  to  A2 . 

D.  Use  CUT  to  generate  an  inequality  which  cuts  off  the  last 
cover.  If  a cuts  have  been  generated  since  the  last  use  of  PRIMAL  1 or  2 , 
go  to  Al;  otherwise  go  to  A2 . 


•'X 


We  call  an  Iteration  of  Algorithm  2 a sequence  of  steps,  including 

{ 

C2,  which  results  either  in  generating  a cut  or  in  fixing  some  variables. 


Numerical  Example 


Consider  the  set  covering  problem  whose  cost  vector  c and  coefficient 
matrix  A are  shown  in  Tableau  3,  and  which  is  obtained  from  the  32-variables 
example  of  [8]  by  removing  9 columns  (12,13,17,20,21,27,28,29,30) 
dominated  by  sums  of  other  columns. 


1 2 3 4 5 6 7 8 9 10  11  12  13  14  15  16  17  18  19  20  21  22  23 
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Iteration  1. 


Al.  PRIMAL  I. 


Step  0.  Rj  « [1 15},  Sx  = 0. 


Step  1.  min  |nJ  is  attained  for  hcT(l)  = {7,13},  N = {7,11,20,22,23}. 
hcR^ 

min  c / l R flM  \ is  attained  for  j(l)  = 23.  S = {23},  R = {5,7,9,12,14}. 

J,Niu) 

Applying  Step  1 four  more  times,  we  determine  the  sequence  of  indices 
j (2)  = 7,  j (3)  = 8,  J (4)  = 16,  j (5)  = 5,  i.e.,  S = {5,7,8,16,23}. 

Step  2.  The  vector  x aefined  by  x^  = 1,  jeS,  x^  = 0,  j et(\S , is  a prime 
cover.  The  associated  upper  bound  is  z^  = 15,  and  T(x)=  M\{l,2,3,8}.  We 


Cl.  STR£?JGTHEN.  There  is  no  cut  yet  to  strengthen. 

TEST  1.  Qq  = {jeN|s^  > 15  - 0}  = 0 (here  s^  = c^ , jcN) ; we  go  to  Bl. 

Bl.  DUAL  1. 

Step  0.  = {l,...,15},  u1  = 0,  s1  «=  c. 

Step  1.  min  _ |N.  ( *=  3,  i(l)  » {7},  u2  = 1 , u2  = u*  i i 7,  and 

ieR1  flT(x)  1 

s2  = (1,1, 1,1, 1,1, 0,1, 2,2,2, 2, 2, 3, 3, 3, 4, 4, 5, 4, 5, 7, 9)  . 

= {7},  M(l)  = ^ = {7},  R2  - R1\M(1)  = {l,...,o,8 15}. 

Applying  Step  1 seven  more  times  we  choose  the  sequence  of  indices 
1(2)  = 13,  i(3)  = 4,  i(4)  = 6,  i(5)  = 9,  i(6)  = 14,  i(7)  = 10,  i(8)  = 15, 

9 

and  obtain  the  vectors  u = (0 ,0 ,0, 1 ,0, 1 , 1 ,0,1, 1 ,0,0,2 , 3,2), 

s9  - (1,1, 1,0, 1,0, 0,0, 0,1, 0,1, 2, 0,1, 0,4, 0,2, 0,1, 0,2),  after  which  Kg  = 0. 

9 9 

Step  2.  u e = 12  > 0 = z , hence  we  set  z = 12  and  stove  s . 

L L 

C2 , TEL.  2.  Qq  = {jcN|s9  > 15  - 12}  «=  {17}.  We  set  = 0,  N - N\{17}. 

Since  Q flS(x)  = 0,  and  E _ s^  = 3 > z - u^e  « 3,  we  go  to  D. 

j«S(x)  J U 


9 

s . 


D.  CUT.  T (x)  = M\{l,2, 3, 8},  S+  = [5,23],  a = 

Step  0.  Wq  = 0,  K1  - <"5, 23],  = 12,  s1  =■  s. 

Step  1.  vt  = mln{ 15-12 , max{l,2}]  = 2,  K*  = {23],  * {13,19,23], 

M(l)  = M23  = {1,... ,4,6,8,10,11,13,15],  anl 

|N  3\£l3, 19,23} | = _win  |n  n[13, 19,23} | = 2. 
icT(x)  D M(l)  1 

Hence  1(1)  = 13,  and  j(l)  = N^OK*  - {23].  We  have  W (13)  = = 

{11,20],  y2  = 12  + 2 = 14  < Zy  = 15,  and 
s2  - (1,1, 1,0, 1,0, 0,0, 0,1, 0,1, 2, 0,1,0, *,0,2, 0,1, 0,0) 

(here  and  In  the  following,  stars  replace  the  entries  corresponding  to 
variables  that  have  been  fixed).  Further,  * {5,23 }\ {23}  *=  {5}. 

Step  1.  v2  = min{l5-14,  l]  « 1 ; K*  « {5};  P,  * {1,2,3,5,10,12,13,15,19,21), 
M(2)  = {5}.  We  have  1(2)  - {5],  j(2)  =■  [5],  W2  = {11, 20) U {14,22},  and 
y3  - 14  + 1 > Zy  = 15. 

Step  2.  We  add  to  (SC)  the  cut 

X11  + X14  + X20  + X22  - 1 

as  the  16th  inequality,  and  store  the  pairs  (W^  = {11,20],  v^  = 2), 

(W2\W1  = {14,22],  v2  = 1). 

Iteration  2. 

A2 . PRIMAL  2 . 

Step  0.  Rt  = {16],  Sl  = {5,7,8,16,23} 

Step  1.  1(1)  = {16},  N1(1)  - {11,14,20,22},  j(l)  = 11. 

S2  = {5,7,8,11,16,23},  R2  = 0. 

Step  2.  The  vector  x whose  support  Is  S = S2\{5]  = {7,8,11,16,23], 
is  a prime  cover,  with  cx  - 16  and  T(x)  = M\{1 ,2 , 3,8, 13} . 
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B2 ■ DUAL  2. 

Step  0.  m = 16.  Initialize  = (0 ,0 ,0 , 1 ,0, 1 , 1 ,0 , 1 , 1 ,0 ,0,0 , 3 , 2 , 1 , ) , 

s1  - (1,1, 1,0, 1,0, 0,0, 0,1, 1,1, 2, -1,1,0, *,0,2,1, 1,-1, 4),  Kj  - (14,22). 

T 

Step  la.  i(l)  = (14),  we  set  u2  = (0,0 ,0, 1 , 0, 1 , 1 ,0 , 1 , 1 , 0 ,0 ,0 , 2 ,2 , 1)  , 


s2  = (1,1, 1,0, 1,0, 0,0, 0,1, 1,1, 2, 0,1,1, *,1,2, 1,1, 0,4),  K~  = 0. 

Step  la.  X°  = {4,6,7,8,9,14,22},  R2  = (3,11,13). 

Step  1.  Since  R2nT(x)  = (ll),  we  set  i(2)  = 11,  = =»in  = l, 

j6Nll 

u3  = u2,  i * 11,  and  s3  = (1,1, 1,0, 1,0, 0,0, 0,1, 1,0, 1,0, 0,1, *,0,1, 1,1, 0,3). 

Further,  F2  = (12,15,18),  M(2)  = (2,4,8,9,11,12,14,15),  and  R3  = (3,13) 

ft 

In  two  more  iterations  of  Step  1,  we  set  i(3)  = 13,  u^  = 1»  anc^ 

i(4)  = 3,  u3  = 1,  and  obtain  the  vectors  u3  = (0 ,0, 1 , 1 ,0 , 1 , 1 ,0, 1 , 1 , 1 ,0, 1 ,2 , 2 , 1) 

and  s3  = (1,1, 0,0, 1,0, 0,0, 0,0, 0,0, 1,0, 0,0, *,0,0, 0,0, 0,1),  at  which  point  R5  = 0. 

Step  2.  u3e  = 13  > 12  = z , hence  we  set  z = 13  and  store  s3  in  place 

L Aj 

- 9 

of  s . 

C2 . TEST  2.  Q-  = (jeN's3  > 15  - 13}  = 0,  2 s3  = 1 < 15  - 13, 

^ i j cS (x) 

hence  we  go  to  A3. 


A3 . PRIKAL  3 . 

Step  0.  Initialize  SL  = {7,8,11,16,23},  Kj  = (1,2,5,13),  ^ = (5,7,9,12,14, 

Step  1.  max  |m,  fl  R.  | « 1 , j (1,1)  - 5;  j (1,2)  = 11.  S = (5,7,8,16,23), 
jcKx  J 

K2  - {1,2,13},  R^,  - (7,9,12,14,16). 

Step  1.  U M OR  - 0;  we  set  R = U M,  = {16}. 

J«K2  J jcS2  J 

Step  2.  -j  (2)  = (20},  S3  = (5,7,8,16,20,23},  R3  = 0. 

Step  3.  U M = M,  S = (5,8,16,20,23).  Since  2 _ s3  = 2 > 15  - 13, 
j cS\(7 } J J*S(x)  J 


we  go  to  D. 
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D.  CUT  generates  the  cut 


Xb  + X11  + X1A  + X19  + X21  + X22 


(which  becomes  the  17th  inequality),  and  the  pairs  (W^  - [b,19,2l),  v^  * 1)  , 
(W2\WL  =*  [11,14,22),  v2  ■=  1). 

Since  a = 2 and  2 cuts  have  been  generated,  we  go  to  Al. 

Iteration  3. 

Al.  PRIMAL  1.  produces  the  prime  cover  x defined  by 
S(x)  = [4,6,7,8,10,11,15,16),  with  T(x)  = M\[3,12,17)  and  cx  = 14.  Since 
14  < = 15,  we  set  z^  = 15,  a 2a  = 4,  and  go  to  Cl. 

Cl.  STFtENGTHEN . Since  z^  was  improved  byl5-14=i,andv  = “ 1 

for  both  cuts  16  and  17,  each  of  them  can  be  strengthened  by  replacing 
W2  with  W . This  yields 

X11  + X20  ^ 1 aad  X6  + X19  + X21  - 1 
as  the  strengthened  inequalities  16  and  17. 


TEST 


_L.  Qq  = [j eNl 3j  > 14  - 13)  = [1,2,5,13,23).  We  set  = x2  = x5  = 


X13  “ x23  = N * [3,4,6 , . . . ,12 , 14, 15 , 16, 18 , . . . , 22) , and  since  x is  still 


a cover,  we  go  to  Bl. 

Bl.  TUAL  1 Droduces  the  vectors  u 

10 


10 


(0,2, 0,1, 0,1, 1,1, 0,1, 0,1, 2, 0,2, 0,0) 


and  s = (*,*,1,0, *,0,0, 0,0, 1,0, 2, *,3,0,0, *,0,0, 0,1,1,*),  with  ue  = 12. 

10 

>ja  - 12  > a wo  Q*r  v-  r=  ^ 

14 

N = [3,4,6,. ..,11,15,16,18,.. .,22).  Since  QQnS(x)  = 0 and 


C2.  TEST  2.  Qq  = [jeNjsj0  > 14  - 12}  = [12, 14}.  We  set  x12  = xix  » 0, 


10 


10 


£ at  = 1 < 14  - 12,  but  E = 4 > 2,  we  go  to  A3. 


j«S(x)  J 


jeN 


A3.  PRIMAL  3.  S (x)  = [4,8,11,15,16,21,22),  cx  = 23. 

D.  CUT. 

Iteration  4 uses  the  sequence  of  steps  A2,  B2 , C2 , D.  Step  C2  sets 


X6  + X19  + X20  ^ 1‘ 


X21  = D generates  the  cut  x18  + X19  + X20  - 1’ 

Iteration  5 consists  of  A2,  B2,  C2,  D,  and  produces  the  cut  X8  + X19  - 1- 


ff  ' * «• 
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Iceration  6.  A2 , Bl,  C2 , x2 > = 0. 

Iteration  7.  A2,  Bl,  C2,  xQ  = xle  = 0. 

Iteration  8.  A2 , Bl.  DUAL  1 produces  the  vector 

u = (0,3,0, 1,2,1 , 1 ,0, 1 ,2 ,3,0 ,0,0,0, 0,0,0, 1 ,0) , with  ue  * 15;  we  set  z = 15, 

L 

and  since  z > z = 14,  the  cover  associated  with  z is  optimal.  This  is  x 
L U (J 

such  that  x j = 1 , j = 4,6,7,8,10,11,15,16,  x^  = 0 otherwise. 

10.  Preliminary  Computational  Experience 
A version  of  Algorithm  1 was  implemented  by  Rohet  Tolani  in  FORTRAN  for  th 
IBM  360/67  computer  at  CMU,  and  a number  of  test  problems  generated  by  Salkir.  a 
Koncal  [9]  were  run.  The  results  are  summarized  in  Tableau  i,  and  compared 
to  an  involutory  basis  cut  algorithm  due  to  Bowman  and  Starr.  The  latter 
is  an  improvement  over  the  Bellmore-Ratliff  procedure,  in  t'.'.at  it  uses  a 
vector  partial  ordering  to  select  a strong  member  of  the  family  of  involutory 
basis  cuts  that  can  be  generated  at  each  step. 


Tableau  1 


No. 

Constraints 

Variables 

Com  ;arison: 

in\  lutory  basis  1 

Cut£ 

1 

104 

133 

18 

15 

120 

2 

200 

300 

35 

30 

49 

3 

200 

413 

55 

49 

365 

4 

200 

500 

238* 

237* 

3000* 

5 

50 

450 

25 

17 

464 

6 

36 

455 

2 

1 

32 

7 

46 

683 

10 

5 

108 

8 

50 

905 

19 

17 

241 

♦Unfinished  run 
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The  discrepancy  between  the  nur-iei  f uer  cion:  and  th-  nu;:.ber  of  cuts 
in  due  to  the  act  that,  if  somic  variable  which  is  in  the  currant  covur  1h 
set  to  0,  the  algorithm  icturns  to  stop  A to  gen*  ate  u net.  cover  instead 
of  generating  a cut. 

As  can  do  seen  from  Tableau  l,  the  algorithm  seems  to  converge  In  a 
surprisingly  smal  1 nut  i itions  (and  cuts;,  though  problem  4 

was  still  unsolved  after  238  iterations,  at  which  point  the  bounds  were 
z = 646  and  zT  = 635  (after  che  first  iteration  the  upper  and  lower  bound; 
were  743  and  597  respectively). 

One  of  ..,e  surprises  >roduced  »j  t test  . 
number  of  variables  that  could  be  fixed,  often  even  in  the 
of  the  tuns.  This  end  the  pattern  of  successive  bound-improvements  is 
illustrated  by  two  typical  r ms  sh  and  ...  ' he  numbers 

ir.  column  one  refer  to  iteiutions  where  either  the  coper  bo„nd,  or  che 
lower  bound,  or  the  number  of  variable  has  changed.  Tb-  u..ueriine  means 
optimum  value. 

Taoleau  2 

froblem  l (104  x 133) 


) Iteration 

zu 

2l 

Variables  1 
left 

1 

1690 

1637 

105 

2 

1663 

87 

3 

1686 

8 

1682 

13 

1678 

14 

1667 

76 

17 

1677 

54 

18 

1 

L_ 

1717 

a. 

i 
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• 1 

f 

I 

\ * 


Tableau  3 

Problem  3 (200  x 413) 


[ Iteration 

Zl) 

z 

Variables 

L 

left 

1 

827 

697 

413 

5 

785 

702 

391 

( 

6 

759 

7 

753 

703 

274 

19 

237 

21 

m 

200 

27 

749 

41 

730 

161 

42 

731 

151 

45 

748 

46 

740 

119 

53 

743 

108 

54 

747 

55 

743 

The  author  is  currently  involved  in  a joint  effort  with  Andrew  Ho 
to  implement  and  test  sever-1  versions  of  Algorithm  2.  Ap_  from  the 
merits  of  various  heuristics,  an  important  aspect  o£  the  ’ reject  concern 


the  relative  merits  of  using  bit  maps  versus  list  structure.  The  result1, 
of  this  project  will  be  reported  separately. 
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*'^A  conditional  lower  bound  on  the  minimand  of  an  integer  progra  s a t .r  ?r 

which  vould  be  a valid  lower  bound  If  the  constrain1  set  we-e  nde 

certain  inequalities,  also  called  conditional.  If  such  : > >r  ' over 

1 bound  exceeds  some  inown  upoer  sound,  then  every  solution  bett  than  e one 

corresponding  to  ;he  upper  bound  violates  at  least  ore  of  the  diti  >r.a  1 

inequalities.  This  yields  a valid  disjunction,  which  can  be  u- 
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the  set  covering  type.  The  family  of  valid  inequalities  derive-  rom 
conditional  bounds  subsumes  as  a special  case  the  Belimore-Ratli: 
inequalities  generated  via  involucory  bases,  but  is  richer  than  the  latter 
class  and  contains  considerably  stronger  members,  where  strength  is 
measured  by  the  number  of  positive  coefficients.  The  paper  discusses  two 
algorithms  based  on  cutting  planes  from  conditional  bounds.  None  of  them 
i uses  the  simplex  method  (though  a variant  based  on  the  latter  is  also 
feasible).  Some  computational  experience  is  presented. 
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